Dynamical behavior of chaos, bifurcation analysis and soliton solutions to a Konno-Onno model

The fractional coupled Konno-Onno model, which is frequently used in numerous fields of scientific and engineering disciplines, is being investigated in the current study in order to gain an understanding of complex phenomena and systems. The two main goals of this study are to be accomplished. Firstly, the research aims to identify novel solitons for the fractional coupled Konno-Onno model using the unified technique, which is currently absent from the literature. Secondly, a novel strategy that hasn’t been previously investigated is phase portrait analysis for both perturbed and non-perturbed dynamical systems. The current study uses appropriate parametric values in phase plane analysis, 2D, 3D, and density plots to ensure the results are physically compatible. The results validate the claim that the technique used in this research to produce complete and uniform responses is not only simple to use and effective, but also substantially faster in computing. The technique is useful for resolving more complex phenomena that arise in engineering and mathematical physics.


Introduction
A potent mathematical tool for modeling and comprehending complicated events characterized by memory effects, long-range interactions, and anomalous diffusion is the fractional partial differential equation (FPDE).These phenomena frequently display behaviors that conventional integer-order partial differential equations are unable to completely explain.The ability of FPDEs to include fractional derivatives, which capture the non-local and non-Markovian nature of the underlying processes, is what makes them special.FPDEs may successfully simulate systems with long-range interactions, where the influence of far-off locations or past events has a large impact on the present behavior, by introducing fractional derivatives [1][2][3].
They are an invaluable tool for simulating and comprehending complicated phenomena in a variety of scientific and engineering disciplines due to their adaptability and capacity to capture nuanced dynamics.Analytical solutions offer precise mathematical formulas that explain how the system will behave in relation to the relevant variables and parameters.The reliance of these solutions on the beginning conditions, boundary conditions, and system characteristics may be learned from these solutions.They provide a more in-depth investigation of the dynamics of the system by revealing the existence of singularities, critical spots, and stability characteristics of the system [4,5].Although numerical methods can be useful in dealing with complicated and nonlinear systems, they may not be sufficient to provide a thorough understanding of the system's behavior.The discretization, numerical techniques, and computational parameters used in numerical solutions generally provide approximations of the results.Analytical solutions, on the other hand, offer a more thorough and accurate description of the system, giving researchers a better understanding of the phenomenon they are studying [6,7].The preference for analytical solutions in the study of FPDEs indicates the need for a deeper comprehension of the properties, dynamics, and underlying mechanisms of the system.Researchers can generate theoretical predictions, gain behavior insightful knowledge of the behavior of the system and direct future research in the area of fractional partial differential equations by locating clear mathematical expressions [8,9].
The coupled Konno-Oono model (CKOS), is an integrable dispersion-less system, that was created by Konno and Oono [21].Both CKOM's integrability traits and its applicability to certain physical events can be credited with its significance.The system's integrability property indicates that its quantities and symmetries are conserved.Depending on the particular system being modeled, this knowledge may have been used in optical fiber and plasma physics.The CKOM is shown as follows in [22]: R xt ðx; tÞ À 2Rðx; tÞMðx; tÞ ¼ 0; M t ðx; tÞ þ 2Rðx; tÞR x ðx; tÞ ¼ 0: The fractional form of the CKOM, i.e., the fractional coupled Konno-Onno model (FCKOM), is preferred in this study because of its advantages in mathematical modeling and optimization.Fractional derivatives provide a more accurate representation of complex physical processes and more accurately reflect system behavior.Fractional derivatives also take into consideration memory effects, in real-world systems and allow the system to retain knowledge from the past.Fractional models are more exact and precise and better suit experimental data than integer-order models.Additionally, utilizing a fractional form enables the application of cutting-edge optimization and control techniques for enhanced system performance.The following is a presentation of the mathematical form of the FCKOM [23]: xt Rðx; tÞ À 2Rðx; tÞMðx; tÞ ¼ 0; at which 0 < δ, γ � 1.The motions of both particles communicating in a medium are represented by the functions Rðx; tÞ and Mðx; tÞ, which have fractional derivatives.Many scholars have previously addressed the CKOM and FCKOM using various analytical techniques; for example, Kocak et al. obtained the traveling wave solution for the CKOM using the modified exp-function technique [24].The tanh-function and expanded tanh-function techniques have been used to create precise outcomes to the CKOM in [25].
The study aims to obtain a novel soliton to an FCKOM via a unified technique (UT) [26] is a novel analytical technique for solving FPDEs.Also, investigate the phase portrait analysis for the perturbed and nonperturbed dynamical system which is obtained by the nonlinear ordinary differential equation (NODE) of FCKOM.This study explores the fascinating area of chaotic behavior in nonlinear systems with the goal of identifying, describing and comprehending the resulting dynamics.The uniqueness of the research, sophisticated analytical methods, and multidisciplinary approach ensure significant contributions to fundamental science and useful applications in numerous fields.Through this work, we uncover new complexities that are at the core of nonlinear dynamical systems, enhancing our understanding of chaotic occurrences and their broader consequences.It applies variable transformations to convert FPDEs into NODEs by using conformable derivatives and then assumes a series-form solution to convert these NODEs into a system of algebraic equations.The families of soliton (also known as solitary wave) solutions for FPDEs are then found by solving the given system of algebraic equations.A soliton is a self-sustaining wave that maintains its shape and speed while neither dissipating nor spreading out.It interacts with other solitons while preserving its distinctiveness and demonstrating stability.Because nonlinearity and dispersion are carefully balanced, solitons exhibit unique properties.To better understand wave behavior, nonlinear dynamics, and integrable systems, they are studied in a variety of fields.They are useful in biological modeling, waterways, and communication networks.Solitons' solutions give insight into how systems work and into the underlying physical processes.
The whole of the work is divided into the following sections.Section, includes the description of the conformable derivative.Section, includes the methodology of the UT.Section, includes the application of the used technique to governing model.The phase portrait analysis is provided in Section.Section, covers the results and discussions.Section, includes the conclusion and future work.

Conformable derivative
The conformable derivative provides an alternative approach to fractional differentiation, which uses the Caputo and Riemann-Liouville procedures [27,28].The conformable derivative has been used in physics, engineering, and data analysis, among other fields.It has been used to represent anomalous diffusion processes, process signals, solve fractional differential equations, and analyze time series.It is significant to remember that the conformable derivative is a relatively recent addition to fractional calculus and that research into its characteristics and applications is still ongoing.A new perspective on fractional differentiation is offered, opening up new possibilities for research and understanding complex processes and systems that display fractional dynamics. Definition: Theorem: Assume p, q : (0, 1) !Q and δ be a differentiable function, then chain rule holds [29] W d ðp � qÞðtÞ ¼ t 1À d qðtÞ dÀ 2 q 0 ðtÞW d ðpðtÞÞj t¼qðtÞ : ð4Þ

Methodology
In order to solve a nonlinear FPDE, the unified technique is used.The UT simplifies the process of finding solutions by providing a consistent approach to addressing multiple types of equations, making it an invaluable tool for researchers in diverse scientific and engineering domains.In numerous instances, the approach has been effectively used to solve fractional nonlinear partial differential equations in physics and mathematics.Let's assume the general form of nonlinear FPDEs: ) is complex valued function and x is spatial and t is temporal term.
Step 1. Assume the transformation of waves are: such as s, q are constant and δ, γ are fractional operator.
Putting the Eq (6) into the Eq (4), the resulting is a nonlinear ODE as: TðJðCÞ; J 0 ðCÞ; J 00 ðCÞ; . ..Þ ¼ 0: Step 2. The solution of Eq (7), can be supposed as following: where l 0 , l m , h m , (m = 1, 2, 3, . .., M) are constants that are (l 0 , l m , h m ), should not be identically zero), M is a term that is to be balanced, by the balancing between the terms that are nonlinear including the greatest derivative and the greatest power of Eq (7) and the function K(C) satisfy the given below differential equation: Where σ is a constant.The resulting solution of Eq (9) are given below that are depends on the σ parameter: (a): When σ < 0, we get hyperbolic solutions K 1 ðCÞ ¼ ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi (b): When σ > 0, we get trigonometric solutions K 5 ðCÞ ¼ ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi (c): When σ = 0, we get rational solution Step 3. Put Eq (8) and second order require derivatives of Eq (8) into Eq (7), now consider the Eq (9), the polynomial which is obtained is power of K(C).
Step 4. Collect all the coefficients of the K(C) which have the same power and furthermore equate each coefficient to zero, then the algebraic system of the equation was derived for l 0 , l m , h m (m = 1, 2, 3, . .., M) and σ.
Step 5. Lastly, solve the algebraic systems of equation by using Wolfram Mathematica and found the parameter values.Putting these parameter values to Eq (8), we get the solution of Eq (2).

Implementation
Let's assume Eq (2) to create an accurate solution.Putting the transformation Eq (6) of a traveling wave into Eq (2), with the help of Eq (4), to get the coupled NODEs are Integrating Eq (11), with respect to C, get Inserting Eq (12), into Eq (10), we get the NODE is where P(C) is complex-valued function, s, q, L are all constants.By using the principle of homogeneous balance, from Eq (13) we derive that M = 1, by a balance between the terms of P 3 (C) and P 00 (C).So that put M = 1 in Eq (8) and solution of Eq (8) is supposed in the form of Derivate Eq (14), twice with the help of considering Eq (9), to get a polynomial in the power of K(C).Collect all the coefficients which have the same powers of the K(C) and after that, set each coefficient to zero.The resultant algebraic system of equation for l 0 , l 1 , h 1 and q as follows When solving the algebraic system of equations, we obtained the solution sets that are Set-1: Set-2: ( By using Set-1 of Eq (16), the solutions for Eq (2) are as follows R 1;1 ðx; tÞ ¼ À ffi ffi ffi L p � ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi sðÀ h 2 À l 2 Þ p À l ffi ffi ffi ffi ffi ffi ffi À s p cosh ð2 ffi ffi ffi ffi ffi ffi ffi À s p ðWÞÞÞ ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi sðÀ ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi By using Set-2 of Eq (17), the solutions for Eq (2) are as follows R 2;1 ðx; tÞ ¼ i ffi ffi ffi L p � ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi sðÀ ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi 2 p ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi sðÀ ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi � ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi ffi  1.

Phase portrait analysis
A mathematical method called phase portrait analysis is used to examine the behavior of dynamic systems.Plotting a series of trajectories on the phase plane, each of which stands for a potential state of the system [30].A dynamical system's behavior can be very dependent on its initial characteristics and conditions, which can result in a variety of phase depictions and trajectory types.Dynamical systems are capable of experiencing events like bifurcation and chaos.Understanding the underlying dynamics and causes of complex systems through the analysis of bifurcations and chaotic behavior aids in deciphering their complex behaviors [31].Following is the phase portrait analysis of the dynamical system for perturbed and nonperturbed systems.

Bifurcation
There are numerous kinds of system bifurcations that can occur.On the dynamics and stability of the system, respectively.Use bifurcation diagrams to see how distinct dynamic regimes transition from one into another.the areas of parameter space where bifurcations take place and the effects they have on how the system behaves.Finding out if there are any higher-order bifurcations will give you a fuller understanding of the system's complexity.

Chaos
An intriguing phenomenon known as chaos can be found in some nonlinear dynamical systems.Even though the system's basic equations are deterministic and clearly specified, it is characterized by sensitive dependency on initial conditions, which causes unexpected and apparently random behavior over time.Numerous scientific fields, including physics, mathematics, biology, engineering, and economics, have conducted substantial research on chaos.
Here are a few of chaos' main characteristics.

Nonperturbed dynamical system
Understanding the qualitative variations in such systems' behavior a parameter or set of parameters is changed is the main objective of bifurcation theory.A nonperturbed dynamical system can exhibit qualitative variations in behavior, and the fundamental mechanisms regulating these changes can be revealed by applying bifurcation theory to the system [32].This analysis contributes to a deeper comprehension of the system's behavior and prospective applications by demonstrating how the system switches between various dynamical regimes as the parameter(s) change.Using bifurcation theory [33], we shall analyze Eq (2) in this section.It is possible to examine Eq (13) as a planar dynamical system by applying a Galilean transformation.
where f 1 ¼ 1 ðsqÞ 2 and f 2 ¼ L ðsqÞ 2 such that s, q and L are parameters.We will now talk about the bifurcations of the phase profiles of system Eq (54), which occur in the parameter space f 1 and f 2 .These are the outcomes of our qualitative analysis.First of all, take note that system Eq (54) has three equilibrium points, which are as follows: Moreover, the system's Jacobian is: Hence, (U, 0) is a saddle for J(U, V)) < 0, (U, 0) is a center for J(U, V) > 0 and a (U, 0) is a cuspidal for J(U, V) = 0.
The following results are potential changes in the parameter values involved.
•  An effective method for comprehending the complex behavior of dynamical systems under parameter variation is through numerical simulation of bifurcation diagrams.Researchers can examine the system's response to changes in the regulating parameter and spot important patterns and transitions as shown in Fig 1 .This approach is frequently used to understand the behavior of complex systems in a variety of disciplines, including physics, biology, economics, and engineering.

Perturbed dynamical system
A dynamical system is said to be perturbed when external disturbances or changes in its properties occur.When perturbations are present, the behavior of the system can be dramatically altered, deviating from that of the corresponding nonperturbed system.Studying how these perturbations affect the dynamics of a system is part of the analysis of perturbed dynamical systems [34].The existence of chaotic behavior can have significant consequences when studying a perturbed dynamical system.The term "chaotic behavior" describes the intricate, erratic, and delicate quality of a system's trajectories, which frequently result from the interaction of nonlinear dynamics and outside shocks.Here is an example of how chaotic behavior might impact how a disrupted dynamical system is analyzed.It is crucial to take into account the unique properties of the perturbations and the system's reaction to them when researching a chaotic dynamical system.To analyze and forecast the system's behavior under perturbations, sophisticated mathematical techniques like nonlinear dynamics and chaos theory are frequently used.A dynamical system's Lyapunov exponents are essential in assessing whether chaos exists there.They offer a numerical assessment of the system's sensitivity to the starting conditions and can show whether the system behaves chaotically.The system shows sensitivity to initial conditions if at least one Lyapunov exponent is positive.An essential trait of chaotic behavior is this sensitivity.Positive Lyapunov exponents reveal the existence of exponential trajectory divergence and the sensitive dependency on beginning conditions, indicating chaos and the system's unpredictable nature.The perturbed dynamical system for Eq (13), as follows: where frequency ψ and amplitude ε.In this investigation, the impact of ψ and ε on the

Results and discussions
It is crucial to assess the improvements and contributions made by the new research by comparing it to earlier work.We can determine the novelty, significance and possible In this research, we obtained novel solitons by using the unified technique.Also, obtained the phase portrait analysis for a perturbed and nonperturbed dynamical system.By employing the UT to FCKOM, get the novel soliton solutions which are hyperbolic, rational and trigonometric.The innovative soliton solutions to the FCKOM that were discovered utilizing the UT may provide important new understandings into a number of physical processes.Understanding wave dynamics, signal propagation, and other nonlinear systems frequently rely heavily on soliton solutions.The novel approaches could lead to the discovery of novel patterns and behaviors in complex systems.By effectively, precisely and adaptably resolving the FCKOM complexity, the UT represents a substantial achievement in the discipline.It is a strong and useful tool for researchers looking at complicated systems including fractional calculus and coupled dynamics because of its capacity to find new soliton solutions, provide analytical insights and handle a wide range of nonlinear equations.Eqs ( 18)-( 25) and ( 36)-( 43) are hyperbolic solutions.Eqs ( 26)-( 33) and ( 44)-( 51) are trigonometric solutions.Eqs (34), (35), ( 52) and (53) are rational solutions.
A specific kind of solution that appears in a number of nonlinear wave equations and demonstrates both soliton and hyperbolic properties is called a hyperbolic solution.Shockwave propagation, conduction of heat, and dispersion problems have solutions that include hyperbolic functions.They offer crucial insights into how these dynamic processes behave, such as the distribution of temperature along a rod.Trigonometric functions are essential for understanding system behavior and designing controls for mechanical oscillations, alternating currents, and voltages.Filtering and frequency analysis are examples of rational functions used in signal processing.They give engineers the ability to work with signals and retrieve crucial information.Applications for these solutions can be found in the study of nonlinear dynamics as well as in the domains of optics, fluid dynamics, and plasma physics.Trigonometric solutions offer an effective foundation for delving into and comprehending a range of mathematical and physical issues.They are a vital instrument in mathematics, physics, engineering, and other disciplines due to their periodic character, association with simple harmonic motion, and applications in numerous fields.For solving equations, modeling relationships, and analyzing systems, rational solutions provide a flexible framework.Their depiction as polynomial ratios allows for a deep understanding of mathematical processes and makes it easier to use them practically in a variety of disciplines, including engineering, physics, statistics and data analysis.

Conclusions
In order to investigate the FCKOM, this study used the unified technique.This method was used to obtain novel soliton solutions, including trigonometric, rational, and hyperbolic ones.Using appropriate parameter values, 2D, density and 3D graphs were used to visualize the results.Additionally, the FCKOM equation's perturbed and unperturbed dynamical systems both performed a phase illustration analysis.Bifurcation theory and chaos theory were https://doi.org/10.1371/journal.pone.0291197.g019also used to analyze the system's behavior.The results of the research showed that the system was composed of both stable and unstable manifolds, bifurcation points and the onset of chaotic behavior.The observed soliton solutions and the analysis of the system dynamics establish a framework for future study in soliton dynamics, stability analysis and the investigation of nonlinear phenomena in related systems.They offer insightful information on the

Fig 1 .
Fig 1. Bifurcation diagram of varying the parameter value of α as shown in Table1.

Figs 15 -
18, represent the phase analysis for an unperturbed dynamical system.Figs 19-22, represent the phase analysis for a perturbed dynamical system.Fig 23, represents the Lyapunov exponent for the existence of chaos dynamical system.

Table 1 . Represent the numerical data of varying α parameter along bifurcation and Lyapunov exponent. TIME a 1 a 2 a 3 Lyapunov exponents Bifurcation analysis
R 2;6 ðx; tÞ